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Abstract 

Many-body techniques for the calculation of quasielastic nuclear matter re- 
sponse functions in the fully antisymmetrized random phase approximation 
on a Hartree-Fock basis are discussed in detail. The methods presented here 
allow for an accurate evaluation of the response functions with little numerical 
effort. Formulae are given for a generic non-relativistic potential parameter- 
ized in terms of meson exchanges; on the other hand, relativistic kinematical 
effects have been accounted for. 
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I. INTRODUCTION 



Quasielastic electron scattering on nuclei has been in the past years the subject of intense 
experimental an d theoretical (see, e. g., Refs. p|-|T7|1) investigations. The first aim of 
the theoretical studies is to test the available nuclear models; once the nuclear physics issues 
are well understood, one might hope to gain insight into other aspects of the problem, for 
instance by extracting with sufficient precision the nucleon form factors. 

In principle, the quasifree regime makes one confident that the physical quantities of 
interest may be computed in a reliable way; in practice, also in this case one has to cope 
with considerable computational problems. Many diverse techniques have been employed in 
the literature. Each of them has its own relative merits and deficiencies and, in general, it 
would be highly desirable to be able to reach some degree of convergence in their outcomes. 

In the following, we shall be concerned with Green's function techniques, as introduced, 
e. g., in Ref. These methods can be, and have been, applied both to finite nuclei 



and nuclear matter, the choice being generally driven by the specific reaction and by the 
momentum regime of interest. Here, we shall focus on the nuclear matter, having in mind 
applications to electron scattering (that is, without the complications introduced by the 
reaction mechanism of hadronic probes) in a range from a few hundreds to several hundreds 
MeV/c of transferred momenta (where the quasielastic peak is sufficiently far from low- 
energy resonances and not too much affected by finite size effects). The use of nuclear matter 
reduces the computational load, thus allowing a more straightforward implementation of 
more sophisticated theoretical schemes: This makes easier to develop and test approximation 
methods that could then be utilized also for calculations in finite nuclei. 

Let us now briefly browse the theoretical framework that we shall discuss in detail in the 
following sections. 

A first choice one has to do in setting up the formalism concerns the treatment of 
relativistic effects. Trivial kinematical effects can be obviously rather important and can 
be included in a straightforward way. The treatment of dynamical effects is more delicate. 
Two main paths have been followed in the literature: Either using field theoretical methods 
(as done, e. g., in the Walecka model and its derivations |0|) or potential techniques (using, 
i. e., phenomenological potentials truncated at some order in the non- relativistic expansion). 
Here, we shall put ourselves on the second path, but, to contain the amount of material, we 
shall employ strictly non-relativistic potentials. The extensions necessary to include higher 
order relativistic terms will be discussed elsewhere (see, however, Refs. [13,^0,21| for a few 
applications). 

Next, one should choose the phenomenological input potential and, in connection with 
this choice, possibly the way of dealing with short-range correlations. All the formulae 
we are going to give in the following sections are based on a generic one-boson-exchange 
potential. They can thus be used both with a bare phenomenological interaction, — such 
as one of the Bonn potential variants, — or with a one-boson-exchange parameterization of 
a G-matrix generated from some potential. The use of an effective interaction derived from 
a G-matrix is a common way of including short-range correlations. One should be aware of 
possible problems due to the use of a local potential to fit non-local matrix elements. At 
least in a few cases discussed in the literature this does not appear to be a reason of concern 



22,23]. On the other hand, possible effects due to the specific quasielastic regime remain 
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completely unexplored: Indeed, G-matrices employed in quasielastic calculations are usually 
generated using bound state boundary conditions, which make them real and practically 
energy independent, while, in general, they might be complex and energy dependent. 

Once we have fixed the effective interaction, we can proceed to consider a hierarchy of 
approximation schemes. 

The lowest order approximation is, of course, given by the free Fermi gas. Then, one may 
include mean field correlations at the Hartree-Fock (HF) level (or Brueckner-Hartree-Fock 
(BHF) if short-range correlations are accounted for). In nuclear matter a HF calculation 
can be done exactly without too many efforts. Nevertheless, we show how a quite accurate 
analytic approximation can be derived, since we shall need the method later to combine the 
HF and the random phase approximation (RPA) schemes. The latter is the last resummation 
technique to be discussed. It should be noticed that even in nuclear matter the calculation of 
the antisymmetrized'RYA response functions is not trivial. Indeed, most calculations that are 
labeled "RPA" in the literature are actually performed in the so-called "ring approximation" , 
where only the direct contributions are kept: In this case, in nuclear matter one gets an 
algebraic equation. Here, we use the continued fraction (CF) technique to provide a semi- 
analytical estimate of the full RPA response (see Refs. || and |J for alternative methods). 
Calculations with this method have been performed both in finite nuclei and in nuclear 
matter p^j25|JT3|j20[| , always truncating the CF expansion at first order, because of the 



difficulty of the numerical calculations. We have pushed the analytical calculation far enough 
to allow not only a fast and accurate estimate of the first order CF expansion, but also of 
the second order one. Since in the CF technique there is no general way of estimating the 
convergence of the series, this is the only way of getting a quantitative hold on the quality of 
the approximation. As noticed before, HF (and kinematical relativistic) effects can then be 
incorporated in the RPA calculation, yielding as the final approximation scheme a HF-RPA 
(or BHF-RPA) response function. 

Of course, many diverse many-body contributions have been left out. It should however 
be noted that the classes of many-body diagrams discussed here, on the one hand already 
allow one to study many interesting features of the quasielastic response; on the other 
hand, the fact of having developed semi-analytical methods reduces to a minimum the 
computational efforts, thus making this formalism a good starting point for the study of 
other many-body effects. 

The paper is organized as follows. In Section 2 the theoretical machinery is set up, 
discussing in separate subsections the treatment of relativistic kinematics and the free, HF 
and RPA responses. The intent of this paper is just to provide theoretical tools, so we do not 
attempt any discussion of the phenomenology of quasielastic scattering. On the other hand, 
in Section 3 calculations based on the formalism previously developed are shown, in order 
to test the convergence of the CF expansion and the importance of antisymmetrization. 
Finally, in the last Section we present a few concluding remarks. 



II. RESPONSE FUNCTIONS 

Let us consider an infinite system of (possibly) interacting nucleons, at some density 
corresponding to a Fermi momentum hp. For the kinetic energies of the nucleons we can 
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choose either the relativistic or non-relativistic expressions, whereas we assume that the 
interactions take place through a non-relativistic potential. For the latter we take the 
following general form in momentum space 

V(k) = V (k) + V T (k)Ti ■ T 2 + V a (k)<n ■ <7 2 + Varik)^ • <7 2 T X ■ T 2 

+V t (k)S 12 (k) + V tT {k)Si2{k) Tl ■ r 2 , (2.1) 

where Sn is the standard tensor operator and V a (k) represents the momentum space poten- 
tial in channel a. Here, we assume that V a (k) has the general form of a static one-boson- 
exchange potential, so that in each spin-isospin channel it is given as a sum of contributions 
from different mesons, V a = X^K^- m the central channels (0, r, cr, err) the contribution 
from any meson can be expressed as the combination of a short-range ("5") piece and a 
longer range ( "momentum dependent" ) pieceQ: 

whereas in the tensor channels (t, tr) is given by 

t4SW^S^_(|^)', *=0,1,2. (2.20 

In Eqs. (|2.2j ), ^\ gy D and cy^N are t ne (dimensional) coupling constant of the z-th meson, 
rrti is its mass and A, the cut-off; to be more general, we have allowed for a choice among 
potentials without form factors or with monopole or dipole form factors. 

Our starting point |26|-p8] is given by the Galitskii-Migdal integral equation for the 
particle-hole (ph) four-point Green's functionF], 



G P aPrts( K + Q,K;P + Q,P) = -G aj {P + Q) G S/3 {P) (2n) 4 S(K - P) 

+iG aX {K + Q)G x , p {K) J ^ri 3 AW (ir + g,ir ; T + g,T)G^, 75 (T + g,T;P + g,p), 

(2.3) 

which is diagrammatically illustrated in Fig. |I[ In ( |2.3| ), G represents the exact one-body 
Green's function, whereas T 13 is the irreducible vertex function in the ph channel. 



The nomenclature stems from the fact that, in the absence of form factors, Vg is a constant and 
is represented by a Dirac <5-function in coordinate space, whereas Vmd is, indeed, the momentum 
dependent piece. 

2 Capital letters refer to four-vectors; small case letters to three- vectors; the Greek letters o,/3, ... 
refer to a set of spin-isospin quantum numbers. 
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FIG. 1. Diagrammatic representation of the Galitskii-Migdal integral equation for the ph 
Green's function, G ph ; T 13 is the irreducible vertex function in the ph channel; the heavy lines 
represent the exact one-body Green's functions. 



Given G ph one can then define the polarization propagat 



or 



= 1 1 i^wr As{K + g ' K] P+Q,n {2A) 

whose diagrammatic representation is displayed in Fig. |2|. Note that for H(q,u) one cannot, 
in general, write down an integral (or algebraic) equation. 

In the case of electron scattering, one can define charge, - - or longitudinal, - - and 
magnetic, — or transverse, — polarization propagators: 

Ul(q 1 u) = tT[6 I I fl(q 1 u)6[] (2.5a) 
U^u) = ^A^U^uj), U^cu) = tr[d^n( 9 ,w)dy (2.5b) 

ij 

where, for brevity, the dependence upon the spin-isospin indices has been represented in 
matrix form, introducing hats where appropriate. In (|2.5|) , I labels the isospin channel and 
the longitudinal and transverse vertex operators are given as follows: 

f Ot° = 1/2 f = at/2 

I Of = r 3 /2 1 = ^ 

The quantity of interest here is the imaginary part of Yl^^^q, uj), since the inelastic scattering 
cross section, — where the momentum q and the energy uo have been transferred to the 
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nucleus, - - is a linear combination of Imil L and ImIT T . It is then customary to define 
longitudinal and transverse response functions 

R L , T (q, = R[^(q, + R[7t(<1, (2-7) 
which are related to IIl,t by 

= -^flf(q,uj)lmUl T ( q ,u J ), (2.8) 

where V is the volume, A the mass number and -p contain the squared electromagnetic 
form factors of the nucleon. 



A. Non-relativistic vs relativistic kinematics 

The response functions introduced above have been defined as functions of the momen- 
tum transfer q and of the energy transfer uo. Actually, it is possible, — and convenient, 
- to define a scaling variable ip, which combines q and uo: This variable is such that the 
free responses in the non-Pauli-blocked region (q > 2kp) can be expressed in terms of the 
unique variable ip (apart from g-dependent multiplicative factors). We shall see that even 
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in the Pauli-blocked region and for an interacting system it is convenient to use the pair of 
variables (q,if)) instead of (q,uj). 

Besides the obvious advantages related to the use of a scaling variable, there is another 
good reason for expressing the responses in terms of ip: In fact, in this way one can define 
response functions that are independent of the form chosen for the nucleon kinetic energy. 
To be more specific: Starting from either a non-relativistic or a relativistic Fermi gas, one 
is always lead to the same expressions for the responses in terms of ((?,V00; to be different 
in the two cases is, of course, the definition of ip in terms of (q,u). 

We shall see in the following subsections that the energy denominators of the free nucleon 
propagators, appearing in the Feynman diagrams for the response functions, can always be 
written as uo — ej^L + , where is the kinetic energy of a nucleon of momentum k and 
k < hp. In the non-relativistic case, one finds 

(0)nr . (0)nr #F / , * & \ , Q n , 

where 

1 / ujm N q\ 

^ = k;[—-2) (2 - 10) 

is the standard scaling variable of the non-relativistic Fermi gas and tun the nucleon mass. 
In the relativistic case, in Ref. |2~9l it had been shown that at the pole it is a very good 



approximation to use Eq. ( |2.9| ) substituting ip nr with 

cum N (l + uj/2m N ) q 



kp 



q 2 



(2.11) 



and multiplying the free response by the Jacobian of the transformation, 1 + uj/itln. 

However, in the calculation of higher order (RPA) contributions, also the real part of the 
energy denominators comes into play and one has to check the quality of the approximation 
far from the pole. With some algebra, — and assuming k 2 /m 2 N « 1, — one can write 



(0)r (0)r 



qkp ip r — q ■ k/kp 



k+q ~ c fc — _ / 
771 TV 1 / 



1 



L) 

1 + 



V 



m N 

qkp ip r — q ■ k/k F 
m N l + uj/m N 



q 2 + 2q-k 



m 



N 



[2.12) 



where, in the last passage, we have replaced the square root with its value at the pole. In 
Fig. ^|, we display the real part of the free polarization propagator (defined in the following 
subsection) using the exact relativistic dispersion relation and the prescription of Eq. (|2.12 ) 



at q = 500 MeV/c and 1 GeV/c as a function of io. We note that the agreement between 
the two ways of calculating Rell( ) is quite good at both momenta. 



Strictly speaking, the validity of this statement is approximate, but quantitatively accurate. 
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FIG. 3. Ren(°) at q = 500 MeV/c and 1000 MeV/c as a function of the energy transfer: Using 
the exact relativistic kinetic energies (solid) and the approximation discussed in the text (dash); 
kp = 195 MeV/c. 



Eq. ( [2.12|) provides an approximation for the free ph propagator: A prescription to get 
the (kinematically) relativistic polarization propagators at any order in the RPA expansion 
(see Section [II D| ) can easily be obtained by noting that U^ n \ — the n-th order contribution 
to the RPA chain, — contains n + 1 ph propagators; then, one has 

Il» r (g, u) = 1 + n^ nr (g, u(l + uo/2m N )). (2.13a) 

V m N J 

Actually, all the response functions derived below are expressed in terms of a generic scaling 
variable ip, as Il^^g, ip): One can then get the non-relativistic response by using the (exact) 
expression ( |2.10|) for ip and the relativistic response by using the (approximate) form (|2.11|) 
and multiplying each polarization propagator by the appropriate power of 1 + uj/misr, i- e. 

1 + nW nr (g,V r ). (2.13b) 

m N J 

Note that in the calculations of Refs. 0,^(J on ly an overall Jacobian factor, 1 +00/171^, has 
been applied to the RPA response functions. 



B. Free response 

Although the free Fermi gas response function is a subject for textbooks (see, e. 



Ref. [@), it is useful to derive it here using a slightly different approach, since it illustrates 



at the simplest level the method we have adopted to overcome a major technical difficulty 
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one meets in nuclear matter calculations, - - namely the presence of 6 functions, which 
considerably complicate analytic integrations. As a side effect, also the calculation of n^°) 
comes out much more compact than in standard treatments. 
From Eqs. ( |2.5|) and (|2.6| ), one immediately finds that 



n 



(0) 

L:I=0 



n 



(0) 

L:I=1 



n 



(0) 

T;7=0 



n 



(0) 

T;7=l 



n<°>, 



where, following ( |2.3|) and ( |2.4|) we have defined 



having set 

Gg(*, 9! v) = -ij ^G®{k + q,k + u)G(°\k, ko) 
G^(k, ko) being the free one-body propagator 

G (0) (Mo) 



9(k - k F ) 6(k F - k) 



h ^ (0) 



7 (0) 
17] K — — 17] 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



The integration over ko in Q2.16Q is straightforward, yielding 



Gg(fc >9 ;«) 



0(fc F - fc)0(|Jfe + q| - fc F ) 0(fc - fc F )0(fc F - |fc + <?|) 



,(o) +M + i71 



+ 



, (0) (0) , ■ 



(2.18) 



which, inserted back into ( P-15 ), would give the standard definition of II^. Instead, let us 



rewrite Gp° as 



d(k F - k)6{\k + q\- k F ) 6(k - k F )6{k F - \k + q\) 



UJ 



(o) 



+ 



, (o) (o) . • 

-^ + e fc+ q -4 +^ 



0(£; F - fc)0(jfc F - \k + q\) 9{k F - k)6{k F - \k + q\) 



(o) (o) . ■ 



, (o) (0) 



(2-19) 



having added and subtracted the quantity in the second line, where we have set rj u 
sign(co>)?7. A few algebraic manipulations then yield 



Hence, from (|2.15|) one gets 



r<°>«. „.,«(**-*) -«(**-!* + «l) 



(0) , (0) . • 



(2.20) 



(27T) 

m N fc| 
~(2^ 



u; — e 



(o) 

fc+q 



+ 



+ ^ 



-a; 



e (0) r 



(0) 



(2-21) 
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Note that only one 9 function forcing k below kp is left, Pauli blocking being enforced by 
cancellations between the energy denominators. In ( 2.21|) , we have introduced q = q/kp and 
the adimensional function 



Q (0) W 



which is easily evaluated, yielding 



1 fi 1-y 2 
dy- 



2J-i ip-y + irju 



(2.22) 



ReQ (0) (V>) =^ + i(l-^ 2 )ln 



(2.23a) 



ImQ {0) (i)) = -sign(w)0(l - ^ 2 )^(1 - ^ 



7l 



sign(c^(l - r)-[P W - P 2 W], (2.23b) 



where P n and Q n are Legendre polynomials and Legendre functions of second kind, respec- 
tively. 

C. Hartree-Fock response 



The HF polarization propagator in nuclear matter is obtained by dressing the one-body 
propagators appearing in n^ ) with the first order self-energy E^ 1 ), so that one can follow 
essentially the same derivation of the previous subsection. The spin-isospin matrix elements 
are the same as for the free response, yielding 



TT HF - TT HF - TT HF 
iL L;I=0 — U L;J=1 ~ iL T;I=0 



n HF 

U T;J=1 



n 



HF 



where 



with e 



(i) _ .(o) 



fc 



+ E«(fc). 



U! 



fc+q 



i (1) I 



+ 



-UJ 



k+q 



JO) 



(2.24) 



(2.25) 



Although the evaluation of the HF response is numerically quite straightforward, in 
Ref. [13] an analytic approximation for Imn HF has been worked out, with the aim of using 



it to include HF correlations in RPA calculations. Here, it will be shown that the validity 
of that approximation is more general, not being limited to the HF response, although in 
the latter case one can directly check the good accuracy of the procedure. 

In any Feynman diagram considered here and in the following, the nucleon self-energy 
enters through the ph energy denominators, 



_(l)nr (l)nr 



-k+q 



+ e 



fc 



OJ 



(k + q) 



k 2 



2m 



N 



2m 



N 



-E«(|fc + 9|) +£«(*), 



(2.26) 



where the non-relativistic expression for the nucleon kinetic energy has been used. In ( 2.26|) , 
one can assume to have always k <kp and \k + q\ > kp. Although the latter may not look 
immediately apparent from, e. g., Eq. ( |2.25| ), remember that cancellations between the 
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energy denominators are such to enforce the Pauli principle; the same will also be true for 
the RPA diagrams^. 

Clearly, if E^(fe) were parabolic in the momentum, the inclusion of the self-energy would 
be achieved simply by substituting m^ with an effective mass. For realistic potentials, a 
parabolic fit of the self-energy over the whole range of momenta is not, in general, a good 
approximation. It is a good approximation, on the other hand, to fit separately the particle 
and the hole part of the self-energy, restricting the fit to the range of momenta actually 
involved in the integration. Since in Eq. ( |2.25| ) (but also in the RPA diagrams discussed 
later) one has k integrated from to k F and, furthermore, \k + q\ > kp, one can set 



k 2 

Z {1) =A + B- , 0<k<k F , 



2m 



N 



(2.27) 



= A + B 



k 2 



2m 



N 



max(g — kp, kp) < k < q + kp. 

Inserting this "biparabolic approximation" back into ([2.26 ) and setting e = A — A one gets 

qk F 
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,(l>r , (l)nr 
t k+q t k 



m 



N 



k 

kp 



(2.28) 



where we have neglected a term proportional to k 2 . It can be expected to be small, since 
k < kp and, typically, q > kp] it has, however, to be checked for any given interaction, since 
it involves the biparabolic fit parameters B and B. In Ref. [0] it had been shown to be 
small for the Bonn potential; the same turns out to be true also for the effective interaction 
employed in the next section. 

Eq. ( |2.28| ) is similar to the expression ( |2.9| ) for the free energy denominator, but for the 
substitutions 



m N 

Vv 



m N 



m N 



1 - 








1 




_ m* N _ 


q 


kp 


q 


2_ 


+ X 
l + B ' 


x = 


1 

kp 



(2.29) 



g^jy _ B q\ 

. q V J 

(or uj — > uj + e). 

In Ref. [|T^| relativistic kinematics had been accounted for by applying the transformation 
uj — > uj(1+uj /2m at) discussed above to the previous formulae. The correct approximation can 
be worked out by starting again from the ph propagator and rewriting it as (AS^^fc, q) = 
£«(*;) -£«(|fc + 91)) 



4 It should also be noted that the infinite Fermi gas is better suited for relatively large momenta 
(q > 2kp), where the conditions above are satisfied by definition. 
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fc+q 



+ e 



(l)r 



u + Jk 2 + m% + ASW(fc, g) + J(k + q) 2 + m 



''N 



CO 



2 + 2c<V/c 2 + + 2(u; 



fc 2 



JV. 



\AEW(k,q) + [AEW(fe,g)] 2 -g 2 



1 + lu/ttin + A^/mjv 



i/j* — q - k/k F 



2qk 

(2.30) 



where 



7?) 



A (1 



v ~~ 1 + 5(1 + cj/mjv + A( 1 )/m JV ) 

^ + x[l + ^(l + w/mjv + AW/2mjv)] 



5 



1 + 5(1 + u/m N + AW/m N ) 
2 _ gfcr 



(2.31) 



2m 



iv 



-x- 



with x already defined in ( |2.29| ). In deriving ( |2.30| ), we have assumed k 2 << m%, evaluated 
the numerator at the pole discarding then any angular dependence and, in the denominator, 
retained only terms at most linear in q ■ kjkp. As one can see, beside the transformation 
uj — > uj(1 + c<j/2mjv) there are other relativistic corrections, both to the effective scaling 
variable and to the Jacobian. 

The quality of the approximations introduced above is good, with at most a few per 
cent discrepancy (except on the borders of the response region, where the Fermi gas is 
anyway unrealistic). Thus, we see that either in the non- relativistic or relativistic case, the 
prescription to include HF correlations in a response function is simply to replace ip with 
ip* and m^v with m* N (and to multiply by a normalization factor when employing relativistic 
kinematics (see Eqs. (|2.13|) ). For instance, from Q2.21| ) one gets 



n HF (g,^) = J 



m 



N 



k'p 



q (2tt) 



[2.32) 



with 



^2.33) 



J r = 1 



Id 



A« 



m N m N 



D. Random phase approximation response 

If in Eq. ( p.3|) one substitutes the irreducible vertex function T 13 with the matrix elements 
of the bare potential, one gets the so called random phase approximation to G ph . In terms 
of the polarization propagator ( |2.4j) one would get an infinite sum of diagrams such as those 
of Fig. |. 
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FIG. 4. Diagrammatic representation of the perturbative expansion for the polarization prop- 
agator in random phase approximation. 



We have already noted at the beginning of Section |I| that, while for the two-body Green's 
function G ph one can introduce an integral equation, this is not possible, in general, for the 
polarization propagator. It becomes possible when one approximates the irreducible vertex 
function T 13 with the direct matrix elements of the interaction. In that case, in an infinite 
system, one gets a simple algebraic equation, whose solution, for the polarization propagators 
JjD and the interaction (|2.1| ), is readily found to be 



TT"ng/ \ n^(g,w) , . 

n x (q, u) = -ma; — zr. r, 2-34 

where ITx^ represents the first order direct polarization propagator: 

</ d =o(i)(<^) = rt \q,tu)W 0{T) (q)rt \q,u;), (2.35a) 
n T;/ d =o(i)(9^) = n(°)( g> w)4[K r(<rr) (g) - V t{tT) (q))n(°\q,u). (2.35b) 

The effect of the exchange diagrams is often included through an effective zero-range inter- 
action, calculated by taking the limit q —>■ of the first order exchange contribution and 
rewriting it as an effective first order direct term ||30|| . Exact calculations, however, show 
that extrapolating this approximation to finite transferred momenta is not always reliable 

A more sophisticated approximation scheme is given by the continued fraction expansion 
PT|,|4]J5]J3^1 . At infinite order the CF expansion gives the exact result as the summation of the 
perturbative series, so that it is not easier to calculate: However, when truncated at finite 
order, it reproduces the standard perturbative series at the same order plus an approximation 
for each one of the infinite number of higher order contributions. The trouble here is that 
there is no general method to predict the convergence of the CF expansion, the only reliable 
test being a direct comparison of the results at successive orders. 

On the other hand, we should note that for zero-range forces the first order CF expansion 
already gives the exact (albeit trivial) result, making one hope that the short-range nature 
of the nuclear interactions allows for a fast convergence. Indeed, all available calculations 
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have been performed truncating the CF expansion at first order P,|5| J^J^JT3|j20|| . Here, as 
anticipated, we shall test the convergence up to second order. 

The CF formalism for the polarization propagator is developed in Ref. [H] for the case 
of Tamm-Dancoff correlations and extended in Ref. || to the full RPA. Instead of following 
the rather involved formal derivation given there, we shall briefly sketch a sort of heuristic 
derivation of the CF expansion (which is, of course, only possible "a posteriori", once the 
meaning of the CF series has been understood). 

Let us assume that we want to build a CF-like expansion for the polarization propagator, 
according to the pattern 

n(°) 

II RPA = . (2.36) 

I — A 

D 

l-C 



We have said that the CF approach at n-th order reproduces the perturbative series at the 
same order and then it approximates higher orders. Thus, if we want to approximate at first 
order in CF the exact RPA propagator (for sake of illustration we drop spin-isospin indices), 



oc 



n RPA = ^ n (n) ) ( 237 ) 



n=0 

we can rather naturally write 



■q(ji) ^ -q(O) 



n(°) 



(2.38) 



where IT^ = U^HVU^ + n« ex is the sum of the direct and exchange first order terms, 
- since this is the correct expression for the direct terms. With this approximation the 
summation is trivial, yielding 

CF1 i - nw /n(°) i - 4\/n(°) - iT 1 )^ /n(°) ' y ' 

We could then add in the denominator of the expression above the exact second order term, 
n^ 2 ), having care of subtracting the approximation to it provided by the first order CF 
expansion, [nW] 2 /n(°). Then, we would get 

TT(°) 



ri - ~ i - nw/nw - {n(2)/n(°) - [n«/n(°)] 2 } 



i - 4vn(°) - n( 1 ) cx /n(°) - {n( 2 ) cx /n(°) - [naw/nw] 2 }' 



(2.40) 



From Eq. ( |2.40| ) it is easy to check that the third order term is approximated as II^ 3 ) 



nW^n^/n'W — [Il^Yn^] 2 }. Then, going ahead in a CF-style expansion we would guess 
for the exact RPA propagator the following expression: 
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X 


u i 


C£ 


c x 


/~1<JT 
U X 


u x 


u x 


L; 1 = 


1 


3 


3 


9 








L; 1 = 1 


1 


-1 


3 


-3 








T;I = 


1 


3 


-1 


-3 


-1 


-3 


T;I = 1 


1 


-1 


-1 


1 


-1 


1 



TABLE I. The spin-isospin coefficients C x (see text), in the longitudinal and transverse 
isoscalar and isovector channels, for the interaction (12.11). 



rr(°) 

TjRPA __ 1 1 



, _ n(1) m( o) nP)«/nw - pi^/nwp ' 

' n( 3 ) ex /n(°) + [n^^/n^] 3 - 2[n( 1 ) ex /n(°)][n( 2 ) ex /n(°)] 

1 n (2)cx /n (0) - [n (1)cx /n (0) ] 2 "' 

(2.41) 

This is exactly the expression one would get from the formalism of Refs. 0,|5j if one had the 
patience to work out the expansion up to third order. Note that we did not assume any 
specific scheme (either Tamm-Dancoff or RPA) in this heuristic derivation. 
Thus, following (|2.34j ) we can write 



n ! PA = r 77T 7-7J, (2-42) 



n(o) 



n w« /n(0 ) _ n^ ex /n(°) 
i - i4 1)d /n(°) _ ng )ex /n(°) - — x } — 

where IIx^ has been defined in ( |2.35| ). Then, truncation at ra-th order requires the calcula- 
tion of the exchange contributions up to that order. 



From ( Q[) we can write 

n£} ex (g, w) = tr[^fl(") ex (g, u)6{] = £ C^U^ an (l: (2.43a) 

on 

ng ex (g,u,) = EVrfO^nW^g,,)^,] = E^/^nS^^o;), (2.43b) 



where the indices a, run over all the spin-isospin channels and the spin-isospin factors are 
condensed in the coefficients C^- an = C^ l] C^ 2) ...C { ^ n) (see Table §. We have introduced 
the "elementary" exchange contribution LT^ ex n containing n interaction lines V ai ...V an , 
namely^ 



5 The following formulae are valid for non-tensor interactions; the treatment of the tensor terms 
is slightly more complex and it is given in Appendix [A]. 
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nitU^) = -* n+1 / |§ ■ • ■ ^?^ )(^ 1 )G^(K 1 + Q)T4 1 (fe a - fe 2 ) ■ ■ • 

• ■ -V an (k n - k n+1 )G^{K n+1 )G^(K n+1 + Q) 

= ( - 1} 7 • • • |^(^-)^ -*■)■•■ 

• • -K„(fc„ - fc n+1 )GS(fc n+1 ,qr;^). (2.44) 

With the definition of given in ( |2.20|) one can, — by a suitable change of integration 
variables, — eliminate all the ^-functions that contain angular integration variables, leaving 
a multiple integral with the following general structure: 



rt£l n M = (-if J ^e(k F - kl ) ■ ■ ■ d ^e(k F - kn+1 ) 



x [ — — V ai (kx - k 2 ) ■ ■ ■ V an (k n - k n+ i 



+ £(c^-u;)]. (2-45) 

In ( |2.45|) , J2{ u — > ~ UJ ) stands for the sum of all the terms generated according to the 
following rules: 

i) take all the terms generated by substituting, in the second line of ( |2.45|) , in one energy 
denominator uj — ► — uj\ then, by doing the same substitution in two energy denomina- 
tors; and so on up to making the replacement uj — ► —uj in all the n + 1 denominators; 

ii) every time (uj - e ki+q + e ki + irj^y 1 is replaced with (-uj - e k%+q + e ki - irju)' 1 , then 
replace ki with — ki — q in the potentials. 

The number of integrations can be reduced by noticing that the azimuthal angles are 
contained only in the potential functions V ai . For typical nuclear physics potentials this 
integration can be done analytically, so that it is convenient to introduce a new function 
representing the azimuthal integral of the potential. To this end, let us define new variables: 



k — k'\ = J k 2 + k' — 2kk'[cos8 cos 6' + sin# sin#' cos(ip — tp')] 



k 2 + k' - 2[yy' + Jk 2 - y 2 ^k' 2 - y' 2 cos(y? - ip 1 )} (2.46) 



x + x' — 2\fx\fx' cos(<y9 — if') + (y — y') 2 , 
where y = k cosO and x = k 2 — y 2 . Then, we can introduce 

W a (x, y; x', y') = C ^V a (k - k') = W Q (x', y'; x, y) (2.47) 
Jo 2n 

and rewrite ( |2.45| ) as 
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■i) 7 



m N 



n+l 



2^ 



2n+2 .1 1 ri-y? 

dyi- 
-i 2 Jo 



1 ^-v 



dx-t ■ ■ ■ \ dy n+1 

l 2 Jo 



n + l 



dx 



n+l 



X- 



ip -yi+ in 



W ai (x 1 ,y 1 ;x 2 ,y 2 ) ■ ■ ■ W an ( 



1 



n+l, 



(2.48) 



For n — 1 one has 



nL 1)ex (g, 



kp 



i 1 /-i-jr 



g J (2tt) 4 7-i y 2 7o 
i 



il /-i-y 



,2 



1 



ip -y + ir] c 



ip-y' + ir) u 



m N 



kp 



q J (2tt)< 



Q ( a ] (o, V) - Qi 1} (g, ) + (o, V + 5) - Qi 1} (-?, V> + ^ , 

(2.49) 



where 



and 



QN(g,V0 = 2 /' dy- t dy'W a "(y,y';q) \ 

J-i ip-y + vq u J-i y-y' 



+ q 



WZ(y,y';q) 



1 /-i-s/ 2 1 /-W 2 
dx 



2 Jo 



dx' W a (x, y + q; x' , y'). 



(2.50) 



(2.51) 



Note that in getting to Eq. ( |2.5U| ) use has been made of the Poincare-Bertrand theorem 
33fl . For the potential (|2.2|) can be calculated analytically (see Appendix [FJ), so that 



the calculation of the first order exchange contribution to the polarization propagator is 
reduced to the numerical evaluation of two-dimensional integrals, — for the real part, - 
and of one- dimensional integrals, — for the imaginary part. 
For n = 2 one has 



m N 



kp r 1 . l r^vi 



f dy^ f 1 dxx f dy 2 \ f * dx 2 f dy^- I 3 dx 
J-i 2 Jo J-i 2 Jo J-i 2 Jo 

1 



1 , 1 r^-vl 



\ q J (2tt) 6 

x 1 , . W a {x x ,yx\ x 2 , y 2 )- , W a (x 2 , y 2 ; x 3 , y 3 ) 

ip-yi + iVw 
( m N 



3 k 6 



ip-y 2 + ir] w ~'~ ip-y3 + ir]u; 

Qi 2 i(0, 0; ij>) - Qi 2 i(0, q- 1>) - Q®(q, 0; ^ + Q<£(g, g; ^) 



V 9 / ( 27 6 

-QS(0, 0; V + q) + 2^(0, -g; ^ + g) + Q®,(-q, 0; ^ + g) - Q®,(-q, -g; V + g)] , 

(2.52) 



where 
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£w(<?i,92;?/0 



i l jo 



dxQ a (x,y + qi,ip + q x ) 



1 



ip - y + ii] L 



and 



W' a {x,y-y') 



dy'~, ; 

i ip — y' + ir) L 



-W^x,y;y') 



1 r l -v' 



dx'W a (x,y;x',y' 



(2.53) 

(2.54a) 
(2.54b) 



For the potential ( |2.2| ) W' a can be calculated analytically (see Appendix [FJ) and one is left 
with the numerical integration of ( 2.53|) and ( p. 54a ), so that the calculation of the second 
order exchange contribution to the polarization propagator is effectively reduced to the 
numerical evaluation of at most three-dimensional integrals. 

Going to higher orders implies a numerical two-dimensional integration for each supple- 
mental interaction line, since, for a potential of the form fl2.2|) , only the azimuthal integration 
can be performed analytically for the interaction lines that are not close to the external ver- 
tices. 

The Hartree-Fock dressing of the nucleon propagators can again be done as explained 



in Subsection [11 Q with the replacements ifj — > ip* and itin 



m 



where ip* 



and m* N have 



been defined in ( (2.29|) and ( p.31| ), multiplying by the correct power of the normalization 
factor 1 + ijjjm^ + /m^ when the relativistic kinematics is employed (see Eqs. Q2.13 ) 
and (12T33|) . 



III. RESULTS 

First of all we have to choose the Fermi momentum. Of course, one could easily perform 
a local density calculation to achieve a better description of finite nuclei: Here, for sake of 
illustration, we prefer to use the pure Fermi gas. The choice of kp can be done in several 
ways: We shall choose an average value according to the formula 

k F = — J drk F (r)p(r), (3.1) 

where p(r) is the empirical Fermi density distribution normalized to the number of nucleons 
and kp(r) = [(3n/2)p(r)] 1 ^ 3 . For 12 C one gets k F ~ 195 MeV/c and this is the value used 
in the calculations that follow. 

Let us start by testing the convergence of the CF expansion. For this purpose, we 
compare in Fig. [5] the longitudinal RPA responses at first and second order in the CF 
expansion using a model one-boson-exchange interaction, V a (k) = <T\ ■ <T2g[in 2 /(m 2 + k 2 )] 
(the spin operators having the purpose of killing the direct (ring) contribution). For values 
of the coupling constant g and of the boson mass m typical of realistic nucleon-nucleon 
potentials one finds that the first and second order results match at the level of a few per 
cent (in the left and middle panels of Fig. §, the solid and dashed curved are actually 
indistinguishable). One has to go to very low boson masses (a few MeV) and, consequently, 
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FIG. 5. Fermi gas longitudinal responses for hp = 195 MeV/c at q = 300 MeV/c, with a 
spin-spin one-boson-exchange interaction, for various values of the coupling constant and of the 
boson mass: Free response (dot), RPA with the first order CF expansion (dash) and RPA with the 
second order CF expansion (solid). Note that in the left and middle panels the dashed and solid 
lines are not distinguishable. The kinematics is non-relativistic. 
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FIG. 6. Modulus of the longitudinal polarization propagator for kp = 195 MeV/c at q = 300 
MeV/c, with a spin-spin one-boson-exchange interaction, for various values of the coupling constant 
and of the boson mass: First order, IlW (dot); exact second order, II^ 2 ) (dash); CF approximation 

^2 



to the second order, n( 2 ) a PP r = IlW /itf ) (solid). 
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to very high values of g in order to find some discrepancies. To understand better these 
results, we display in Fig. |6] the modulus of the polarization propagator at first order, 
(dot), at second order, iT 2 ) (dash) and the approximation to n^ 2 ) generated by the first 
order CF expansion (see Section pDf), n( 2 ) appr = /U^ (solid). From inspection of 
the curves, it is clear that the first important element to guarantee a good convergence 
is the range of the interaction: Indeed, for m = 800 MeV (short-range) II^ 2 ^ and n( 2 ) appr 
practically coincide independently of the strength of the interaction. This, of course, should 
be expected, since for zero-range interactions the first order CF expansion gives the exact 
result. For masses of the order of the pion mass one starts finding discrepancies between 
iT 2 ) and ]T 2 ) appr : However, for realistic values of the interaction strength the second order 
contribution turns out to be one order of magnitude smaller than the first order one, thus 
making these discrepancies having no effect on the full response functions (Fig. ||). 

To understand these results it may be useful to compare the strength of the interactions 
employed here to the one of one-pion-exchange, g^m 2 = / 2 /3 = 0.33 (in natural units). 
With the same units, the cases with m = 100 MeV correspond to gm 2 = 0.26 and 0.65; 
those with m = 800 MeV to gm 2 = 16.7 and 41.7; for m = 1 and 10 MeV one has gm 2 = 1.3 
and 0.65, respectively. 

To summarize, from the left and middle panels of Fig. |6] one can understand that the 
validity of the CF expansion originates out of the interplay between range and strength of the 
interaction: For short-range potentials, — where the conventional perturbative expansion 
may not converge, - - the CF technique gives a good approximation of the propagators 
at all orders; for long-range (on the nuclear scale) forces, the CF approximation is less 
accurate, but the relative weakness of the interaction already guarantees the convergence of 
the conventional perturbative expansion. One has to go to unreasonably low masses to find 
a situation where the interaction range is very long and and are of the same order 
(right panels in Fig. |6]). 

We can thus conclude that accurate calculations of nuclear response functions in the 
antisymmetrized RPA can be performed at first order in the CF expansion. The same 
conclusion is supported also by calculations with a realistic effective interaction, — such 
as the G-matrix parameterization discussed below, - - and including HF and relativistic 
kinematical effects. 

Finally, it is interesting and important to test the validity of the ring approximation, - 
where exchange diagrams are not included, — since this approximation has been widely used 
in the literature because of its simplicity. In this scheme, the effect of antisymmetrization 
is mimicked by adding to the direct interaction matrix elements an effective exchange con- 



tribution (see, e. g., Ref. |30|]). For details see also Ref. |34|, where a prescription, suitable 
for the quasifree region, to determine the effective exchange momentum has been given. 

In Fig. [7|, then, we display the ring and RPA responses of 12 C at q = 500 MeV/ c, using the 
G-matrix parameterization. It is apparent that the only channel where the ring approxima- 
tion works reasonably well is the spin-isovector one, — which, incidentally, is the dominant 
one in (e,e') magnetic scattering; it is less accurate in all the other channels, particularly 
in the scalar-isoscalar one. The same considerations apply also when HF correlations are 
included in the ring and RPA responses. Note that these results confirm those of Ref. || , 
where a comparison of ring and RPA calculations had been done using a numerically rather 
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FIG. 7. Fermi gas responses for fcp = 195 MeV/c at q = 500 MeV/c, with the G-matrix 
parameterization discussed in the text: Free response (dot), ring approximation (dash) and RPA 
(solid). The kinematics is relativistic. 



involved finite nucleus formalism, 
been employed. 



Also in that calculation the G-matrix of Ref. 221 had 



IV. CONCLUDING REMARKS 



In this paper we have illustrated a fast and compact scheme for the calculation of the 
fully antisymmetrized RPA response functions in nuclear matter, based on the CF expan- 
sion. The fast convergence of the CF series for typical nucleon-nucleon potentials has been 
demonstrated, thus making this technique a very convenient tool for the exact resummation 
of the RPA diagrams. On the other hand the poor performance in most spin-isospin channels 
of the ring approximation, — where the exchange diagrams are not included, — has been 
confirmed. Accurate approximations for the inclusion of the relativistic kinematics and of 
HF effects have also been discussed and tested. 

Although other classes of contributions are also necessary to make contact with the 
electron scattering phenomenology, — such as meson exchange currents and higher order ph 
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configurations, — we believe that the methods discussed in this paper provide a good, - 
because of the accuracy, — and convenient, — because of the simplicity, — starting point. 
As mentioned in the Introduction, it would also be interesting a comparison, under the 
same approximation schemes, with other approaches, such as those based on the relativistic 
models of nuclear structure. 



APPENDIX A: TENSOR INTERACTION IN THE EXCHANGE DIAGRAMS 

The n-th order exchange polarization propagator in presence of tensor interactions has 
an expression slightly more complicated than (|2.48 ), because the tensor operators do not 



allow, in general, for a factorization of the azimuthal integrations. A generic diagram with 
m non-tensor and n — m tensor interaction lines can instead be written as 



ra+l / , \ 2n+2 



n<-£L,. ..,,(?.") = (-D" (^) (|) 

x / dyi\ [ 1 dxi ■ ■ ■ [ dy ri+ i\ [ " +1 dx 
J-i 2 Jo J-i 2 Jo 



•Emi Vrm x m+l j Vm+X j 

1 

x W / a m+ l...a„(Wl)2/m+li •••) x n+l, Vn+l)' 



( A1 ) 



where W ai has been defined in Q2.47| ) for the non-tensor channels and 



' 27r dip m+1 f 27T dip n+ i 



2tt 



Wa m+1 ...a n ( x m+l, Um+l] %n+l, Vn+l) = 2™ m } J > , A?'» / — o "■ / 

. . , i Jo 2tt Jo 

x ^«m+i(^m+i ~~ k m+ 2)Sii 1 {k m+ i — k m+ 2) ■ ■ -V an {k n — k n+ i)Si n _ m j{k n — k n+ i). (A2) 



In the last expression we have introduced the tensors 

Sij(k) = 3kikj - Sij, (A3) 

such that J2ij <Ti(TjSij(k) = Suits). 

The first order case is rather simple, since one gets again ( |2.49|) -( p.51 ) with 



W a ix,y;x',y') = [ * -^V a {k - tf)S„{k - k'). (A4) 

JO 27T 

At second order, however, one can use Eqs. ( J2.52p - fl2.54p only when just one tensor interac- 
tion is present. 
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APPENDIX B: FIRST AND SECOND ORDER EXCHANGE DIAGRAMS 



We give here the explicit expressions for the first and second order exchange diagrams, 
based on the potential fl2.1|) - (pl2|) . In ( |2.49| ) and Q2.50| ) we have seen that 



ni 1)ex (g, 



where 



to 





\ 2 k% 




1 (2tt) 4 



Qi 1} (o, v») - Qi l> (g, V) + Q ( a'(o, 1> + q) - Q ( *\-q, 1> + q) 



>(i) 



dy'W a "(y,y'; q) 



y-y' + q 



:bi) 



(B2) 



whereas from ( 2.52 ) and Q2.53 ) one has 



m N 



U6 



Q<£,(0, 0; V) - Q^(0, g; V) - Q^(g, 0; V) + Q^(?, g; 



,(2) 



(2) ^ 



)( 2 ) ^ ~ , 



V ? 7 W 

-Qi 2 i(0, 0; ^ + g) + 2^(0, -g; ^ + q) + Sg^-g, 0; ^ + g) - Q^-g, -g; ^ + g)] , 

(B3) 

where 

(2) _ f 1 1 Z" 1 ^ 2 1 

Qaa'iQ^q^;^) = dy- dxg a (x,y + g x ;^ + gi) , _ . g a '{x,y + q 2 ;^ + q 2 ) 

(B4) 

and 



G a (x,y;ij}) = / dj/' 



1 



i *p -y' + ir]i 



-W^x,y;y'). 



(B5) 



For a meson-exchange potential the quantities that can be calculated analytically are those 
given by Eqs. (EE37D , (gg), (|2T54bD and (|2T5l| ), namely 



Jo 2n 

r 27T d(f 



and 



W' a {x,y-y') 
W^(y,y'; q) 



o 2vr 



(non-tensor) 
y Q (fc - k')S zz (k^k') (tensor) 



1 r 1 -?/ 

2 Jo 

1 /-1-2/ 2 1 /-l-'i/' 2 



2 Jo 



drfW tx {x,y',rf,'if) 

dx'W a (x,y + q; x',y'). 



dx - , 
2 Jo 



(B6a) 
(B6b) 

(B7) 
(B8) 



In any channel a the potential is expressed as a combination of the terms displayed in fl2.2|) . 
Then, for each of them one finds 
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W s (x,y;x',y') 



W M D(x,y;x',y') 



W T N(x,y;x',y') 




1 

2 

e = o 



(B9a) 



(B9b) 



95, i = 

g s (X 2 - /i 2 )w a (X\x,y;x',y'), £=l 
g s (\ 2 - fi 2 fw b (\\x,y;x',y'), £ = 2 

' 9md ti 2 w a (ti\x,y;x',i/), 
g MD /j, 2 [w a (iJ,\x,y,x',y') - w a (\\x,y;x',y')], 
9md fi 2 [w a (n\x, y; x\ y') - w a {\\x, y; x\ y') 

-(\ 2 -H 2 )w b (\\x,y;x>,y% £ = 

' 9TN{[3(y ~ y') 2 + /J> 2 }w a (/j,\x,y;x',y') - 1}, 
fi-TN{[3(y - y'f + fi 2 ]w a (/j,\x, y; x', y') 
-[3(y-y') 2 + \ 2 }w a (\\x,y;x',y')}, 

gMiHy - y'f + A 

x [w a (fi\x, y; x', y') - w a (\\x, y; x', y')\ 
-(A 2 - /i 2 ) [3(y - y') 2 + A> b (A|x, y; y')}, I = 2, 

where again £ labels the power of the form factors, we have introduced the adimensional 
form factor cut-off, A = A//cp, and meson mass, \i = m/kF, and we have defined 



(B9c) 



w a (X\x, y- x', y>) = {[A 2 + (y - y'f + x + xf - Axx'}^ 2 

A 2 + (y - y') 2 + x + x' 



w b (\\x,y;x',y 
For W' one finds 



{[X 2 + (y - y') 2 + x + xf - Axx'} 3 / 2 ' 



W' 5 {x,y-y') 



WM D (x,y;y') = < 




1 

2 



W^(x,y;y') 



f^(l-y' 2 )/2, £ 
\ g s (\ 2 - fi 2 )w> a (\\x,y;yi), £ 
{g 5 (\ 2 -Li 2 ) 2 w b (\\x,y;y>), £ 

' g M D^ 2 w' a (^\x,y;y'), 
9md /x 2 K(/i|x, y; y') - w' a (\\x, y; y')], 
9md /i 2 K(/i|x, y; y') - w' a (\\x, y; y') 
-(A 2 -/i 2 K(A|x,y;yO], 
f ^ TN {[3(y - y'f + Aw'M^ y\ y') - (1 - */' 2 )/2}, 

^TN{[3(y-y') 2 + /^ 2 K(/^, y; y') 

-[3(y-y0 2 + AV«(A|x,y;yO}, 
^TN{[3(y-y') 2 + V 2 }[wMx,y;y') - w' a (\\x, y; y')\ 
-(A 2 - /x 2 ) [3(i/ - y'f + A 2 M(A|z, y; y')}, 



£ = 
1 

£ = 2 



(BlOa) 
(BlOb) 

(Blla) 
(Bllb) 



£ = 
£= 1 
€ = 2, 



(Bile) 



where 



<(A|x, y;y') = - 



x In 



A 2 + (y - y'f + l-y' 2 -x+ ^[A 2 + (y - y') 2 + 1 - y' 2 + x} 2 - 4(1 - y' 2 )x 



w' b {\\x,y;y') = - 



A\ 2 + (y-y'f 



2 [A 2 + (y - y'f] 

\ 2 + (y-y'f-l+y' 2 + x 
0A 2 + (y - y') 2 + 1 - y' 2 + a;] 2 - 4(1 - y' 2 )x 



(B12a) 



(B12b) 
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Finally, for one finds 



(B13a) 



W^ D (y,y';q) 



f gs[(l-y 2 )/2}[(l-y' 2 )/2], £ = 
W^y,y';q) = { g s (X 2 - ^)w^X\y,y f ;q), £=l 
[ 9s(X 2 - fi 2 ) 2 w'b(\\y, y'\ q), 1 = 2 

g M nV 2 wZ(iJ,\y,y';q), £ = 

g M -D^ 2 [wMy,y';q) -v%(\\y,y';q)], £=l 
g MDf i 2 [wMy,y'-q)-w';(\\y,y';q) - (X 2 - ^ 2 )w^X\y,y';q)}, £ = 2 

(B13b) 

{ ^T N {[3(y - y' + q) 2 + /i 2 ]<(/i|y, y'; q) 

-l(l-y 2 )/2][(l-y' 2 )/2]}, £ = 

grn{[3(y -y' + qf + fi 2 ]wMy,y'; q) 
-[S(y-y' + q) 2 + X 2 ]w':(X\y,y';q)}, 

gM[3(y -y' + qf + v 2 } K0%, J/'; o) - v'\ <?)] 

{ -(A 2 - fi 2 )[3(y -y' + qf + A 2 ]<(A|i/, y'; q)}, £ = 2, 



W^ N (y,y';q) 



l (B13c) 



where 

w:(X\y,y';q)= 1 -{ 

-[A 2 + (y - y' + qf + 2 - y 2 - y' 2 \ - 2(2 - y 2 - y' 2 ) In |2[A 2 + (y - y' + qf]\ 
+\/[A 2 + (y - y< + q) 2 + 2 - y* - y' 2 ] 2 - 4(1 - ^)(1 - y' 2 ) 
+2(1 -y 2 ) 

A 2 + (y-y' + ?) 2 + z/ 2 -y /2 + ^ 2 + (?/-2/ / + ?) 2 + 2-y 2 -y' 2 ] 2 -4(l-y 2 )(l-y' 2 ) 



x In 



+2(1-/) 



x In 



X 2 + (y-y' + q) 2 -y 2 + y' 2 + J[X* + (y-y> + qf + 2-y*-y> 2 f-±(l-yi)(l-y> 2 ) } 

(B14a) 

A 2 + (y - j/' + qf + 2 - - y' 2 - 



x- 



A 2 + (y - ^ + qf + 2 - y 2 - y' 2 ] 2 - 4(1 - y*)(l - y' 2 ) 
X* + (y-y> + qf 

(B14b) 
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